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1. Introduction 

Thanks to the exact chiral symmetry, lattice simulation with the overlap fermions [|I|] provides 
a unique opportunity to approach low energy hadron physics. One of the major issues in this 
energy region is the consistency between QCD and chiral perturbation theory (ChPT). With the 
overlap fermion, the continuum ChPT can be applied without modification, that is not true with 
other fermion formulations that violate either chiral or flavor symmetry. 

In this talk, we present two-flavor QCD simulation using the overlap fermion and the Iwasaki 
gauge action. On a 16^ x 32 lattice, we generate 10,000 trajectories [^, ^ ^, ^ ^ at six different sea 
quark masses. We explicitly suppress zero-modes of hermitian Wilson-Dirac operator H\f = /jDvi? 
in order to avoid discontinuity of the overlap-Dirac operator Dqv = — 75 • sgn(//vi?(— nio))) 
by introducing extra Wilson fermions which are irrelevant in the continuum [^]. As a result, the 
topological charge of our gauge configurations is fixed. We have chosen the trivial topological 
sector Qtop = 0. 

In Section 2, we present the spectrum calculation focusing on new techniques using eigen- 
modes of the dirac operator. After making corrections to our data due to the finite size effects 
(FSE) in Section 3, we carry out chiral extrapolation of the corrected results using the NNLO 
ChPT formulae and compare the results with phenomenological values. 



2. Spectrum calculation 

In advance of the spectrum calculation, we compute and store 50 pairs of the lowest-lying eign- 
modes of Dqv on each gauge configuration. With these data, we decompose the quark propagator 
as 

50 u'Jx)-u'Jy) 

sMy)ai^= I (2-1) 

where the indices a,j8 represent spin x color and / labels eigenmode. m° is the bare quark mass. 
While we use the conventional CG algorithm to obtain the second term with the source vector 
with low mode contributions projected out, we construct the first term from the eigenvalue A; and 
eigenvectors u'fy,{x). Since it is the low mode contribution which is computationally dominating the 
CG iteration, Sq^^ can be obtained with much lower computational cost ( x 8 less). 

Computation of the meson correlator is done on 500 configurations separated by 20 HMC 



trajectories. When decomposing quark propagator as in (2.1), a part of the meson correlator is 
composed purely from the low mode contribution. By replacing this part by its average over the 
source location (low-mode-averaging |^ ^), we can improve the data. Its advantage is apparent in 
Figure |l] where the effective mass of pseudo-scalar correlation functions with (filled symbols) and 
without (open symbols) low-mode-averaging are compared. 

We obtain pion mass and decay constant from the pseudo-scalar correlators with smeared 
source and with local source. By fitting them simultaneously with single exponential functions 
sharing the same mass, we obtain the matrix element of the local operator | (0 \P\ 7i) p, from which 
fji is computed, with a better statistical quality than solely using the local-local correlator. is 
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Figure 1: Effective mass from the P-P correlation function with smeared source for the lightest three quark 
masses Wsea = Wvai- Open symbols represent data from conventional correlation function while filled ones 
are from low-mode-averaged correlator. 



obtained though the axial Ward-Takahashi identity 

A = 2m°(0|P|7r)/m2, (2.2) 

without further renormalization. Note that we are using the /^j = 130 MeV normalization. 

Throughout this calculation, we use ro = 0.49 fm as an input. Our lattice scale at the chiral 
limit is 



a = 0.1184(12) (11) fm. (2.3) 

By using its central value oT^ = 1.67 GeV, the pion mass data points are in 290 MeV <mji< 750 
MeV. Our heaviest m° is roughly corresponding to the strange quark mass m^. 

We are interested in the chiral extrapolation of m^/niq and where is the quark mass 
renormalized by the renormalization factor 

Z,^(2GeV) = 0.742(12), (2.4) 

which is obtained non-perturbatively through the RI/MOM scheme on the lattice |^. 



3. Finite size corrections 



3.1 Standard FSE 



Some of our data points are in the region mj^L < 3.0 where finite size effect could be significant. 
Therefore, it is important to estimate the correction to our data. To this end, we use the result of 
analytic calculation by Colangelo et al. pO|]. Developing from the Liischer's formula [11], which 



relates the mass-shift in a finite box to the amplitude of scattering of pions, they calculated the 
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finite volume corrections for niji and fj^. Thougii their calculation includes NNNLO of ChPT for 
and NNLO for /;;:, we apply their NNLO results for both quantities. The NNLO effects depend 
on the low energy constants (LECs) Ti, 2,3,4 of the Nf = 2 ChPT. At the scale of the physical pion 
mass m^^^ = 139.6 MeV, they are estimated [12] as 



;-phys ^ _o.4±0.6, (3.1) 

//^y^ = 4.3 ±0.1, (3.2) 

^-phys ^ 2.9±2.4, (3.3) 

[/■^y^ = 4.4 ±0.2. (3.4) 



We use these phenomenological values to correct our data. The errors in ( [3.1[ )-(|3.4|) are reflected 
in the results assuming a gaussian distribution. Because the Liischer's formula is based on the full 
theory, we can use this result only for the data with the same mass for valence and sea quarks. 

3.2 Correction for fixed 2top 

Since our numerical simulation is done at the topological charge fixed to zero, our observables 
are not free from the artifact. However, for large enough volume, local quantities such as hadron 
mass or matrix elements do not depend on Qtop- In fact it is verified by a saddle-point expan- 



sion [ |13| ] that the difference between the correlation function at fixed topological charge and that 
in the 0-vacuum is of ff{V^^). By using this relation to pion mass and decay constant, it can be 
shown that the leading correction is proportional to their second derivative with respect to 6 depen- 
dence at the saddle point ds = i'Gtop/ (VXt)- Xt is the topological susceptibility, which is calculated 
on the same set of configurations [jlj]. At NLO of ChPT we obtain the corrections 



m^^" ^ 1 



m^{d=0) \6Vx, 



(3.5) 



1+77777 ( ) iini i (3.6) 



Md=0) 4Vx, V 471/ J \ \ mf'' 



where m^^^^{d)^ = 2BomqCos{d /Nf) is the tree-level 0-dependent pion mass. The fixed Qtop cor- 
rection starts at the tree-level for nij^ while it does at NLO for fj^. 

In Figure we illustrate the effects of finite size effects for m^/niq (left) and fji (right). In 
each panel, the original data (black diamonds) are corrected by the standard FSE (blue squares) 
and further by the fixed topology effect (red circles). For the data with the lightest quark mass, the 
standard FSE is 4.5 % and 6.0 % for m^/rriq and fj^, respectively. For rn^/mq in small mass region, 
two finite size effects appear with an opposite sign to almost cancel each other 

4. Chiral extrapolation 

(0 \ tree 
m^j = 

IB^niq, we chose t, = {mn/ {^nfn))^ expansion parameter This choice enables us to fit the 
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Figure 2: Results of finite size corrections for m^/m^ (left) and f;^ (right). Data points as a function of ms, 
are slightly shifted for visibiUty. 



ChPT formulae (almost) independently to m| / and /;;: as a function of which is a measurable 
quantity at each quark mass. The NLO formulae are 



{ml/mq)^^^ = 2Bo (1 + § In^ - ., 
(A)nlo = /(l-2^1n§+2/4§) 



(4.1) 
(4.2) 



where Z3 4, Bo and / are to be fitted. 4 are LECs at the scale 47if. Since these two formulae do 
not share any fit parameter, independent fits are possible. 

We also examine the NNLO formulae which are written as 



{ml/mq) 



NNLO 



2Bo 



l+^ln^ + ^(^ln^)2-Z3(^-9^2i„^)^J^;^_4^;phys^l6) j^2i^^ 



(4.3) 



NNLO 



/ 



l-2§ln^+5(^ln^)2 + 2/4(^-10^2ij^^)^3 /~phy,^^|^2ij^^ 



(4.4) 



where coefficients of the term a\ and a2 contain unknown low energy constants associated with 
the NNLO counterterms. In carrying out the chiral fit with these formulae, an input 



;~phys_7;-phys^g^-^phys_^5j^ 



(4.5) 



is introduced from phenomenological estimates (3J) and ( [3.2[ ) with ff^^^ = 130.0 MeV. 

In our target range < 1^ ^ 0.1, an approximation ^^In^ —2.5^^ is numerically rather 
precise. Therefore, it is reasonable to consider yet another fit ansatz with modified NNLO formulae. 
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Figure 3: Results of chiral extrapolation for m^/mq [GeV] (left) and fj^ [GeV] (right) as a function of ^. 
For each panel, three kinds of fit (NLO, NNLO and NNLO') by using all data points are shown. 



which we call NNLO' 
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l + ^ln^ + -(^ln^)2-Z3^ 



+ a;^2^^(^3), (4.6) 



(A)nnlo' =/[l-2^1n^+5(^ln^)' + 2/4^]+a^^2 + ^(^3). (4.7) 

In Figure |3| we show the results of chiral fits to the FSE corrected data of m|/m^ and fjt, 
respectively. In each panel, we compare NLO (red curves), NNLO (blue) and NNLO' (green). 
Note that we cannot compare the quality of the fit with NNLO and others by the values of ;^^/dof 
because we carry out the simultaneous fit for the former and independent fits for the latter. 

Figure ^ compares quantities from our fit: /, £ = Bq • /^/2, l^'^^ and 1^^^ . In the figure we 
also plot phenomenological values of / [ 15 ] and l^^^^ (([3 .31) and (^)) as well as the result of our 



previous calculation in e-regime ||T1: £1/3 =0.251(7)(11) GeV. 

We see that, for l^"^^ , large error of the phenomenological estimate covers all of our results. 
For other quantities, results from NLO fit are inconsistent with the NNLO results. This implies the 
failure of the NLO formulae to describe the data in the quark mass region up to m.,. Results from 
NNLO and NNLO' are consistent with each other and with phenomenological estimates. 

For the future we are plamiing to extend the present analysis to thk and fK on the Nf = 2+1 
dynamical lattices (see [ p^ for the status). Numerical check of the estimate of FSE in a bigger 
volume is also planned. 



Numerical simulations are performed on Hitachi SRllOOO and IBM System Blue Gene Solu- 
tion at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale 
Simulation Program (No. 07-16). This work is supported in part by the Grant-in- Aid of the Ministry 
of Education (Nos. 13135204, 15540251, 17740171, 18034011, 18340075, 18740167, 18840045, 
19540286, 19740121, 19740160). 
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Figure 4: Comparison of physical quantities obtained from chiral fit ansatz: / (upper left), E^/^ (upper 
right), f^^^ (lower left) and F^^^ (lower right). In each panel, red circles and blue squares are corresponding 
to fits with 6 and 5 lightest data points, respectively. 
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